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the scattering kernel for the collision integral truncated after the second term, we derive the diffusion limit 
for the Boltzmann equation by truncation of Chapman-Enskog or Hilbert expansions with reaction and 
Ph ■ collision scaling. We also give asymptotically sharp results obtained by the use of an additional time scaling. 

The diffusion limit determines the diffusion source in the Isotropic Diffusion Source Approximation (IDSA) 
\ of Boltzmann's equation [30l [5] for which the free streaming limit and the reaction limit serve as limiters. 

' Here, we derive the reaction limit as well as the free streaming limit by truncation of Chapman-Enskog 

or Hilbert expansions using reaction and collision scaling as well as time scaling, respectively. Finally, we 
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CN ■ Introduction 



Some astrophysical processes like nuclear burning in the sun, neutron star mergers, or the collapse 
and explosion of a massive star, emit enormeous quantities of neutrinos from captures of electrons 
. or positrons on nucleons. If such events occur in our Galaxy, the neutrinos can be recorded on Earth 



>< 



together with the electromagnetic spectrum of the corresponding astrophysical events. However, 
since the neutrinos only weakly interact with matter, they are more difficult to detect than photons. 
On the other hand, due to their weak interactions, they can escape from much denser and hotter 
matter than the photons so that they carry information about this extreme state of matter directly 
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to the terrestrial observer. Moreover, they are not obscured by dust and may even be detectable 
from events that are not visible in the electromagnetic spectrum. By the detection of neutrinos, 
models describing nuclear burning in the sun could be confirmed, neutrino oscillations were detected 
that lead to extensions of the standard model of particle physics [3l|T8], and the general core-collapse 
scenario of type II supernovae was confirmed by the detection of neutrinos from supernova 1987A, 
i.e., the explosion of a massive star in the close- by Large Magellanic Cloud |22l, !§]. 

Astrophysicists started to appreciate the role of neutrinos in the explosion of massive stars 
in 1966 already [Ml |T], only 10 years after the postulated neutrino has been detected as a free 
particle [24J. Ever since, simulations of core-collapse supernovae attempted to take the drain of 
energy and lepton number by the emission of neutrinos into account [9l [Tj [37t I46j . In fact, more 
than 99% of the released gravitational binding energy of the neutron star that forms at the center 
of the stellar collapse is carried away by neutrinos [6]. It turned out that neutrinos are not only 
important as messengers from density and temperature regimes that are not otherwise accessible by 
observation, but they also have a large impact on the dynamics of the supernova itself. The neutrino 
energy is considered crucial for reviving and feeding the shock that finally leads to the explosion of 
the star. But the exact mechanism of this process is still not completely understood. The main issue 
is that sufficiently accurate and reliable radiative transfer methods were initially only implemented 
in spherical symmetry, as for example the solution of the 0{v/c) Boltzmann equation [321 HOl H3] 
or the complete general relativistic Boltzmann equation [28^ |42] . But the explosion dynamics and 
the interaction of the neutrinos with the convective stellar matter and magnetic fields at the core 
of the explosion are fundamentally three-dimensional processes [21^ I23[. It is therefore crucial to 
investigate the reliability of efficient approximations of the radiative transfer equations that can 
more efficiently be used in three-dimensional models. For instance, one seeks approximations of 
the Boltzmann equation that capture the main processes of neutrino transport in different regimes 
while being computationally more economic. Well-known approximations like the Multi-Group 
Flux Limited Diffusion approximations [101 EH] and the Variable Eddington Factor Method [3H [36] 
have successfully been extended to axi-symmetry. 

The main physical behaviour of neutrinos in core-collapse supernovae is that they are trapped 
by the matter in high density regimes, where their dynamics is governed by diffusion, whereas they 
are practically freely streaming in low density regimes further away from the core. In both cases, 
the Boltzmann equation can be reduced considerably. This is the underlying idea of the decompo- 
sition used in the Isotropic Diffusion Source Approximation (IDS A) of Boltzmann's equation [30] , 
Concretely, the distribution function of the neutrinos is decomposed additively into a trapped and 
a streaming particle component on the whole domain. The resulting equations that govern these 
two particle components are reduced to the main physical properties of the components. The 
major challenge of the approximation is to find an appropriate coupling between the two reduced 
equations. In the IDSA they are supposed to be coupled by a source term. The concrete form of 
this source term, the diffusion source, is based on the diffusion limit and, for non-diffusive regimes, 
limited from above and below on the basis of the free streaming and reaction limits. The derivation 
of the diffusion limit in terms of an asymptotic expansion is sketched in |30t App. A] and given in 
full detail in [5]. 

The main purpose of this paper is to derive the diffusion, free streaming and reaction limits in the 
IDSA by Chapman-Enskog and Hilbert expansions of variations of the reaction and collision scaled 
as well as the time scaled Boltzmann equation. We give the concrete order of the approximations 
to the Boltzmann equation in terms of the scaling parameter e > and thus provide a deeper 



2 



understanding of the approximation quality for radiative neutrino transfer given by the IDSA. Here, 
we concentrate almost entirely on the analysis of the IDSA. Computational aspects, discretization 
and numerical solution techniques for the IDSA and the Boltzmann equation accompanied with 
numerical results which compare the two models are presented in [301 [5] . 

Concretely, the main topics and results discussed in this paper are as follows. In Section [H we 
recall two models used for the simulation of core-collapse supernovae, namely a fully coupled model 
of hydrodynamic and 0{v/c) Boltzmann equations and the Isotropic Diffusion Source Approxima- 
tion. We restrict our considerations to spherical symmetry although extensions to the full 3D case, 
in particular to an IDSA model in 3D, are feasible [301 Sec. 4]. With regard to the IDSA, we explain 
how the hypotheses on the trapped and streaming particle components lead to the corresponding 
reduced equations. We perform a Legendre expansion of the scattering kernel that is crucial for 
the derivation of the diffusion and the reaction limit of Boltzmann's equation. The three regimes 
of the IDSA represented by the diffusion source and its limiters are discussed mathematically and 
physically. Finally, we address some mathematical issues of the IDSA. 

Section [21 which contains the asymptotic analysis, is the core of this paper. Here, we perform a 
Chapman-Enskog expansion of Boltzmann's equation and then consider Hilbert expansions. The 
results can be summarized as follows. Introducing a scaling parameter e > that accounts for 
large reaction and collision terms, we can derive the diffusion limit of the Boltzmann equation up 
to the order 0{e^) if we neglect the same terms that needed to be neglected in the "leading order" 
approximation in |30[ [5]. With the same scaling and without the need to neglect these terms, the 
reaction limit of the Boltzmann equation up to the order 0(e) can be derived. Here, the limit used 
in the IDSA is obtained if the free streaming component decreases strongly enough asymptotically. 

In principle, the same results can be derived by Chapman-Enskog expansions as well as by 
Hilbert expansions. However, the Hilbert expansion approach seems the more natural one since in 
fact, it is finally used for the actual comparison with Boltzmann's equation. For the application 
of Chapman-Enskog expansions, we resort to a lemma that provides approximations of Chapman- 
Enskog expansions by Hilbert expansions. This lemma requires additional assumptions on the 
asymptotic behaviour of components of the expanded distribution function in the Chapman-Enskog 
expansion that are thus not needed if we start with the Hilbert expansion. Nevertheless, we find 
that these connections between Chapman-Enskog and Hilbert expansions give a good insight into 
the nature of these expansions that are interesting on their own. 

As the main sharp asymptotic result, we provide a derivation of the diffusion limit of Boltz- 
mann's equation and thus of the diffusion source up to the order 0{e^) with the use of an addi- 
tional time scaling by which it is no longer necessary to neglect any further terms. Finally, the free 
streaming limit of the Boltzmann equation is derived up to the order 0{£) if only time scaling of 
the equation is considered. The paper closes with a motivation why it is natural to use the free 
streaming and the reaction limits as limiters of the diffusion source as it is done in the IDSA. 

1 Models for Core-Collapse Supernovae 

In this section, we introduce two models for the simulation of core-collapse supernovae. The first 
model consists of a coupled system of hydrodynamic and Boltzmann equations. The second model 
is the Isotropic Diffusion Source Approximation (IDSA) of the first model and has been developed 
in Liebendorfer et al. [30] . Although the two models can be considered both in spherical symmetry 
and without symmetry in 3D, we restrict ourselves to the spherically symmetric case here. The 
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description of the models resembles in parts the presentation given in however, the discussion 
of them is more comprehensive here. 



1.1 Fully coupled model 

A widely accepted model for the simulation of core-collapse supernovae is a coupled system of the 
hydrodynamic equations ([T]) for the background matter (including leptonic number conservation) 
with Poisson's equation ^ for gravity and Boltzmann's equation ([3]) for the radiative transfer of 
neutrinos. We refer to papers of different groups [9l[33], [41] . |43j and [42] who have implemented 
different algorithms for the same physics and, in particular, to the comparison paper j29j . However, 
all these approaches are restricted to spherically symmetric cases. The presentation in this paper 
is based the version used in |30j . We write the system of equations in short as 

G(u,<A) + S(u,/) (1) 
AttGp (2) 
i(u)-x(u)/ + C(u,/) . (3) 

V ' 

Jif) 

This system is given in spherical symmetry, i.e., for {r,t) G fi-r := (0, i?) x (0,T) for some radius 
R> and some end time T > 0. Here, (0, i?) represents the spatial domain U C M'^, which is the 
ball around the origin with radius R, where the matter of the star that needs to be considered is 
thought to be contained. In the following we specify the equations in more detail. 



^ + V • F(u) + Vp(u) 

Idf df ^ , ,df ^ , ,df 

^ V ' 

»(/) 



1.1.1 Hydrodynamics, state equation and gravity 



System ([T]) contains the hydrodynamic equations for an ideal fluid, i.e., the Euler equations [25], as 
well as a balance law for the leptonic number. The state vector u € M^, which is a function on 
the nonlinear flux function F : M'^ — )• and the source terms G = G(u, </>) and S = S(u, /) have 
the form 
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while, as usual, the gradient of the term p(u) = (0,p(u), 0, 0)"^ with the pressure p(u) only enters 
the momentum equation. The variables contained in the state vector u are the density field p 
and the velocity v in radial direction of the background matter. More precisely, since the baryon 
number is the conserved quantity in astrophysical hydrodynamics, whether relativistic or classical, 
p = nbrrit, is considered to be the product of the baryon number density rib and the baryonic 
mass constant rrib- Furthermore, the total energy density E of the matter is given as the sum of its 

2 

internal energy density pe and its kinetic energy density p^. The pressure p is assumed to be given 
by an equation of state p = p{p,'d{e),Ye) where -i? is the temperature (depending on the internal 
specific energy e) and Ye = ^ is the net electron fraction of the matter [301 p. 1179]. Here, Ue is the 
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electron minus the positron number density, so that pY(, = Uemi, is the net electron density upto the 
constant rub [Ml PP- 640/1]. The appearance of Ye in the equation of state (in [30] the Lattimer- 
Swesty equation [26] is used) takes the compositional degrees of freedom into account. Under the 
extreme densities in the stellar core and in case of nuclear statistical equilibrium that we consider, 
the pressure is essentially given by the degeneracy pressure of the electrons. This dependency of 
the pressure on Yg necessitates to add a fourth balance law (leptonic number conservation) to the 
Euler equations which describes the evolution of Yg- 

On the right hand side, the bilinear term G(u, (/)) accounts for the effect of the gravitational 
field of the matter onto itself, see, e.g., [441 Sec. 7]. The Newtonian gravitational potential (p on Q 
is given by the Poisson equation ([2]) where G denotes the gravitational constant [25', p. 6]. With 
regard to equations ([I]) and ([2]) we recall that in spherical symmetry, with the radial unit vector e^, 
the gradient and the divergence operators are given by 



for scalar fields V : K ^ M and vector fields F = (F^, 0, 0)^ : M ^ (see, e.g., [Ml P- 679]). 

The coupling of system ([T|) for the matter with the Boltzmann equation ([3]) for the neutrinos is 
provided by the source term S(u, /) which adds to global momentum, energy and leptonic number 
conservation ([I])2~(IIj4- Generally speaking, S(u, /) consists of sums of certain energy and angular 
moments of the right hand side of ([3|) and is nonlinear in u and, due to the restricted set of neutrino 
interactions with matter in Section [1.1.21 linear in /. For instance, if the neutrinos were isotropic, 
Sv{u, f) would become the negative gradient of the neutrino pressure which corresponds to the 
pressure gradient of the matter on the left hand side of ([l|), compare [301 Eq. (24)]. For further 
details we refer to [33] and [3^ . 

Remark 1.1 With regard to supernova simulations, it is crucial and, for our purpose, enough to 
know that the right hand side S(u, /) of the hydrodynamic equations does not depend explicitly on f 
but only on energy and angular moments of f. In fact, the explosion of a star can be seen by the 
velocity of the background matter and can be inferred by the detection of neutrinos fjT^ . However, 
only the neutrino flux and energy spectrum are observable, but not their angular distribution f . 
It therefore suffices to consider only moments of f instead of f itself. This is done both in the 
IDS A and in the asymptotics that we perform for the radiative transfer equations in Section [H 

Both hydrodynamics ([1]) and gravitation ([2|) could be replaced by more accurate general rela- 
tivistic versions [27^ 141] . For instance, the gravitational binding energy released in a core-collapse 
supernova is equivalent to about 1/10 of the mass of the remaining neutron star. This fact had 
already been claimed by W. Baade and F. Zwicky two years after the discovery of the neutron |2j. 
Recall that the mass conservation equation ([T])i is in fact a baryon number conservation equation, 
so that the mass defect is not ignored here. Since the gravitational binding energy is almost com- 
pletely converted into neutrino energy [Jl], the energy balance equation ([TJs can be regarded as the 
conservation of energy including gravitational binding energy. The supernova models considered 
here are also valid if only a protoneutron star forms at the core and later collapses into a black hole 
due to continued accretion before, during or after the explosion, when part of the material does 
not reach escape velocity and falls back. 




and 
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1.1.2 Radiative transfer 



The radiative transfer equation used in [30] for the transport of neutrinos is the 0{v/c) Boltzmann 
equation ([3|) in spherical symmetry. It represents both the special and the general relativistic 
transport equation for massless fermions moving with the speed of light c up to the order 0{v/c) 
where v is the velocity of the background matter. We refer to |35j Sec. 95] for a derivation and [H] 
for a discussion of possibly additional 0{v/c) terms that might have to be considered. In practice, 
the fraction v/c can be expected to be around 1/10 in most of the computational domain except 
for the shock region where v/c ~ 3/10 can occur and at very high densities deep inside the neutron 
star where it gets very close to 1, see |27j . In that sense the 0{v/c) Boltzmann equation must be 
regarded as an approximative model to the physical reality. 

The Boltzmann equation ([3]) describes the evolution of the neutrino distribution function 

/: [0,T] X (0,i?] X [-1,1] X (0,^] ^ [0,1], (t, r, ^, a;) ^ /(t, r, /i, . (5) 

Here, /i € [—1, 1] is the cosine of the angle between the outward radial direction and the direction 
of neutrino propagation and uj S (0, -E] with some maximal > is the neutrino energy. In the 
formulation of ^ that we use, both phase space variables ^ and uj are given in the frame comoving 
with the background matter, cf. [33*, p. 640]. In more general approaches, one could replace (0, i?] 
by (0,oo). 

With regard to the notation we use the Lagrangian time derivative 

d d d , . 

in the comoving frame. Furthermore, the coefficient function -F^(u) weighting the //-derivative 
in ([3]) can be decomposed as 

where the term Fj^^ accounts for the change in propagation direction due to inward or outward 

movement of the neutrino. The second term F^{u)^ represents the angular aberration (i.e., the 
change observed in the comoving frame) of the neutrino propagation direction due to the motion 
of the fluid. Finally, the product of the coefficient function 



2 j dlnp 3v\ V 
^ cdt cr J cr 



ijj 



with ^ in ([3]) accounts for the (Doppler-) shift in neutrino energy due to the motion of the matter. 
The left hand side of the Boltzmann equation is abbreviated by where P is a linear operator. 

The dependency of T) on the hydrodynamic variables in u occurring in the comoving frame 
vanishes in case of a static background with frozen matter where one can pass to the laboratory 
frame by setting -F^(u) = Fi^{\jl) = 0, i.e., 

c at or r Ofx 
and using Lorentz transformed quantities in ([7|), see [35, §90, §95]. 
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Although in the infall phase [9l p. 787] it is enough to consider only electron neutrinos fe, for 
postbounce simulations [30^ p. 1179] one needs at least two Boltzmann equations in order to obtain 
the transport of both electron neutrinos and electron antineutrinos v^- In general, one needs to 
include muon and tau neutrinos and their antiparticles, too [H]. With our set of interactions given 
below, all different types of neutrinos are transported independently, so that one has to deal with 
up to six Boltzmann equations that are all coupled with the hydrodynamic equations by the source 
term 5(u, /). Since the basic structure of these equations is the same, it is enough for our purpose 
to consider only one Boltzmann equation ([3]) as a prototype. We point out, however, that these 
equations are coupled if we include interactions between different neutrino types like neutrino pair 
reactions [H p. 774/5] or, as a future perspective, neutrino flavour oscillations [3t ITS], 

As for system ([T]), the coupling with hydrodynamics is provided by source terms on the right 
hand side of Their values usually differ for different neutrino flavours and their antiparticles. 
The interaction of neutrinos with matter includes emission and absorption 

+ n ^ p + Ve 

of electron neutrinos v^. or electron antineutrinos by protons p or neutrons n, the forward reactions 
being known as electron e" or positron capture. Analogous reactions occur in case of electron e~ 
or positron capture by nuclei. They depend on the state of the background matter and result in 
neutrino emissivity j (w, u) and absorptivity x(lj, u) whose sum is the neutrino opacity x = J + X in 
([3]), cf. [33l p. 640] and [30l p. 1177]. Formulas for the nonnegative expressions u) and xl*^, u), 
that are nonlinear in uo and in u, for both electron neutrino and its antiparticle have been derived 
in [91 pp. 822-826]. 

By the term C(u, /), the right hand side of ([3]) also accounts for isoenergetic scattering of 
neutrinos (or antineutrinos) on protons, neutrons and nuclei. It is a linear integral operator in /, 
the so-called collision integral, and reads 



C(u, /(t,r,^,Lj)) 



6^2 



c{hc) 



-1 



R{u,n,fi',uj){f{t,r,^i',uj) - f{t,r,n,uj))d^j! 



(9) 



where the isoenergetic scattering kernel i?(u(t, r), //, /i', w) is symmetric in ^ and ^' and depends 
nonlinearly on all its entries, see [91 pp. 806/7, 826-828] for concrete formulas. In Q Planck's 
constant is denoted by /i, the term corresponding to f{t,r,fj,',uj) accounts for in-scattering while 
the term corresponding to f{t,r,n,uj) represents out-scattering [30l p. 698]. Note that if / does 
not depend on we have C{f) = 0. As an immediate consequence of the symmetry of the collision 
kernel with respect to fi and fi' we obtain for any / 

1 

C{u,f)dfi = 0. (10) 

1 

Further neutrino interactions with the background matter as considered in [9l p. 774] such as, 
e.g., neutrino electron scattering are neglected in |30] . In the following we omit the dependency 
on u in the terms of the Boltzmann equation and abbreviate the right hand side of ([3]) by j -|- J7'(/) 
where i7(/) is linear in /. 
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1.2 Isotropic Diffusion Source Approximation (IDSA) 

In this section we give a short introduction to the Isotropic Diffusion Source Approximation (IDSA) 
that has been developed in [30]. The aim of this approximation of Boltzmann's equation ^ is to 
reduce the computational cost for its solution, making use of the fact that ^ is mainly governed 
by diffusion of neutrinos in the inner core and by transport of free streaming neutrinos in the 
outer layers of a star. The following ansatz for the IDSA intends to avoid the solution of the full 
Boltzmann equation in a third domain in between these two regimes as well as the detection of the 
corresponding domain boundaries. 



1.2.1 Ansatz: Decomposition into trapped and streaming neutrinos 

We assume a decomposition of 

f = f + r (11) 

on the whole domain into distribution functions /* and supposed to account for trapped and for 
streaming neutrinos, respectively. 

With this assumption and using the linearity 'D{f) = V{f) + V{f^) as well as J{f) = J{f) + 
J'(/*), solving the Boltzmann equation ([3|), i.e., 

V{f)=3+J{f) 
is equivalent to solving the two equations 

v{f)=j + j{f)-i: (12) 

v{n = j{n+i: (13) 

with an arbitrary coupling term S(r, t, ^, w, u). For the IDSA one establishes approxima- 
tions of these two equations arising from physical properties of trapped and streaming particles, 
respectively, and one determines an appropriate coupling function E(r, ^, w, u). Follow- 
ing Remark II. 1^ the aim is to deal only with angular moments of /* and instead of the full 
distribution functions themselves. 



1.2.2 Hypotheses and their consequences for the trapped and streaming particle 
equations 

Concretely, one uses the following hypotheses. First, the trapped particle component /* as well as 
S are assumed to be isotropic, i.e., independent of /i. Taking the angular mean of equation ()12p 
this leads to the trapped particle equation 

idinp a/* . , 

— r "I ; — — = J ~ XJ ~ ^ ) (14) 

cdt 3 cdt ouj 

in which we slightly abuse the notation 

f = \ j\ f'dn (15) 

concerning the domain of definition of the isotropic /* given in ([5|) . The isotropic source term E is 
treated in the same way. Here, and in what follows, we always assume that we can interchange the 



8 



integral and the differentiation operators. Details of the computation can be found in the proof of 
Lemma 12.61 

Next, /* is assumed to be in the diffusion limit, which is physically at least justified for the 
inner core of the star. In order to derive the diffusion limit, a Legendre expansion of the scattering 
kernel R{fi, fj.' ,uj) with respect to its angular dependence, truncated after the second term, is used 
in [30|. App. A] for an approximation of the collision integral, see Subsection 11.2.31 In fact, this 
approximation is essential for the derivation of the diffusion limit and thus the resulting definition 
of the source term E. In [5], this derivation based on a Chapman-Enskog-like expansion of ([3]) 
is performed in detail. In case of very high densities, even the diffusion term can be neglected. 
Then /* is in the reaction limit and satisfies a reaction equation. It is the purpose of this paper 
to give derivations of the diffusion, the reaction and the free streaming limit, to which we turn in 
the following, by mathematically concise applications of Chapman-Enskog expansions and Hilbert 
expansions, see Section [2j 

Secondly, the streaming particle component is assumed to be in the free streaming limit. 
This justifies to neglect the collision integral in (I13p . which by (llOp . however, vanishes anyway after 
angular integration. Furthermore, it justifies to neglect the dynamics of background matter so that 
one can use the laboratory frame formulation ([7]) of ([3]) with frozen matter as the left hand side 
of (|13p (here, we also neglect the Lorentz transformation). For the same reason one can assume the 
free streaming particle component to be in the stationary state limit and drop the time derivative 
in ([7]) which then, after angular integration, becomes the streaming particle equation 

Due to dH) this equation can be reformulated as a Poisson equation for a spatial potential ip of the 
first angular moment of 

For the solution of ()16p . the approximate relationship that has been suggested by Bruenn in|27). 

between the particle flux and the particle density of streaming neutrinos is used. Here, Ry{uj) > 
is the energy dependent radius of the neutrino scattering spheres. In addition to spherical sym- 
metry, this approximation is based on the assumption that all free streaming particles of a given 
energy uj are emitted isotropically at their corresponding scattering sphere [30.3 P- 1178]. Hence, 
the flux i f_-^ f^{uj)^d^ can be approximated in terms of the particle density | f^{oj)d^ by the 
geometrical relationship ()17p . The geometrical factor that correlates these terms is the only place 
in the IDSA where the angular distribution of neutrinos, namely of the free streaming neutrinos, 
is (approximatively) encoded. This factor is equal to 1/2 if r < R^{uj), which is a consequence 
of the isotropy of / inside the scattering spheres, and increases up to 1 in the limit r — > oo. The 
latter expresses the fact that the neutrinos tend to stream radially outwards so that the distribution 
function / accumulates at /Li = 1. 

1.2.3 Legendre expansion of the scattering kernel 

As mentioned above, we now seek for an approximation of the collision integral suited to derive 
the diffusion limit and also the reaction limit of Boltzmann's equation that we envisage in the next 
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section. The approximation is performed by a Legendre expansion of the scattering kernel. For 
an introduction to Legendre expansions by spherical harmonics we refer to [45^ pp. 302, 391-395]. 

2 

Concretely, the Legendre series for ■^^^^R{fi, fi' ,u}) reads 



— — 3i?(^,M',a;) = — ^(2/ + l)0i(a;) / Pi{cose)dip, (18) 
c(. cj vr Jo 

with the Legendre polynomials P/, Z = 0, 1, . . ., where 9 is the angle between the incoming and the 
outgoing particle and 

cos((9) = + [(1 -//2)(l -^'2^]5cOSV9. (19) 
With the first two Legendre polynomials given by 

Po (cos e) = l, Pi (cos 6*) = cos 6* , 

truncation of the series after the second term gives the approximation 

^2 ^ / 1-211 r2Tr \ g 

■^^^j-^R{fi,fi\u) — i(l)o{uj) 1 d(p + 3cl)i{uj) cos{e)dip] = -(/)o(cj) + -0i(cj)/x/x' . 

Inserting this into the collision integral ([9]), without explicitly mentioning the dependency on t, r 
and oj, one obtains 

C{f) « 1 1'^{4>0 + 30iM/i')(/(^') - /(/u))d/u' = -M + <t>ol fdfi + 3^01^ ffsdfi , (20) 

which is an affine function in /x expressed in terms of /, the zeroth and first angular moments of/ 
and the opacities (po and (pi. Together with the emissivity and absorptivity of neutrinos by the 
matter on the right hand side of ([3]) the latter gives rise to the definition of the neutrino mean free 
path 

A — ^ - ^ (21) 

This definition is motivated by the fact that A/3 occurs as the diffusion parameter in the diffusion 
limit of the Boltzmann equation that will be derived in the following section by asymptotic expan- 
sions (see (j63p ). It is clear that the smaller the diffusion parameter is, the smaller the diffusion 
of neutrinos represented by the diffusion term in (j63p becomes which physically corresponds to a 
smaller mean free path. 

Finally, we remark that a collision kernel which can be expanded as in (jl8p with (jl9p is always 
symmetric in fi and /i' since cos(^) in ([19]) has this property. For the same reason, (pOj) as well as 
any truncation of the Legendre series in p8p is symmetric in /i and fi'. In particular, property (jlOp 
also holds for the approximation of the collision integral in (j20p . 



1.2.4 Diffusion source and system of coupled equations 

Using the truncated Legendre expansion of the collision kernel as in the previous subsection, the 
isotropic diffusion source S = Si^s is derived in [301 App. A] by a Chapman-Enskog-like asymp- 
totic expansion. The derivation that provides the diffusion source in a certain "leading order 
approximation" is explained in detail in [5j. Here, we quote the result which we will reobtain 



10 



by Chapman-Enskog and Hilbert expansions in Section [2] where we will also explain the precise 
meaning of the "leading order approximation" . The diffusion source reads 




We emphasize that, in spite of (jlOp . the influence of the collision integral is contained "in leading 
order" in S (^oj, /*, ^ J^-^ f^dfi, uj by its dependence on u which contains the mean free path A and 
also the emissivity j and the opacity x- 

By physical reasons outlined in [30l p. 1177] and also justified in Subsections 11.2.61 and II .2771 
the coupling term S is limited from above by j and from below by 0, i.e., for the coupling of the 
trapped particle equation (fH|) and the streaming particle equation HfTU^i the source term 

S := min {max [Sids , 0] , j} (23) 

is used in [30]. As a result, one obtains the coupled system of equations (fill) . (fT6]l . (i23|) as the 
Isotropic Diffusion Source Approximation IDS A of Boltzmann's equation ([3]). 

The limiters T, = j in (jl6p and S = in (jl4p can be regarded as representing the free streaming 
limit and the reaction limit, respectively, of the Boltzmann equation ([3]). In Section [21 we will 
also derive these limits by Chapman-Enskog and Hilbert expansions under certain assumptions. 
A mathematical motivation why these limits enter formula (I23p in the form of limiters for the source 
term S will be addressed at the end of this paper (see Lemma l2.23p . 

As a consequence of the limiters j and in (I23p , one can consider three different regimes of the 
system (jl4p , ()16p and (j23p depending on the three alternatives S can assume in formula ()23p . We call 
these regimes according to the corresponding limits of the Boltzmann equation ([3]) and emphasize 
that transient situations that lie in between the limits are treated in the IDSA as belonging to one 
of the three regimes. In this sense the terms "regime" and "limit" have different meanings in the 
following. 

1.2.5 Diffusion regime: S = Sids 

Including the diffusion term in ()22p . noting its formulation in spherical symmetry in (j4]), into the 
left hand side of ([H]) provides a diffusion-reaction-type equation 

for /* balanced by the emissivity j and an effective absorption — x (^/* + ^ J^-^ f^dfij of trapped and 
streaming neutrinos. Consequently, in the diffusion regime there is only absorption but no emission 
of streaming particles. However, trapped particles can be converted into streaming particles. The 
equation for the latter reduces to 

i.e., — f/* is a potential of the streaming particle density. In the diffusion limit, where one can 
assume r < Ry \n pT]) . this flux is just half of the particle density. The diffusion regime turns out 
to be the only regime where the equations for /* and i.e., (p^ and ([25]) . are effectively coupled. 
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Since the basic ansatz (jlip - (jl3p for the IDSA assumes an additive couphng term E in the 
equations for the trapped and the streaming particle component on the whole domain, i.e., also in 
regimes where ^ cannot be supposed to be in the diffusion limit, the diffusion regime of the IDSA 
is likely to be used where the diffusion limit does not apply. For more extreme cases, i.e., where 
the limiters in (j23p apply, the diffusion source is altered in order to account for other limits of the 
Boltzmann equation to which we turn now. 



1.2.6 Free streaming regime: S = j 

If we neglect the dynamics of the background matter, i.e., if we assume p to be constant in (I14p . 
this equation has a solution /* that decreases exponentially if we have S = j and even more so if 
S > j. Therefore, the trapped particle component /* vanishes quickly in the case S > so that 
its equation (|14p as part of the Boltzmann equation ([3]) can eventually be neglected. However, 
neglecting the right hand side of (|14p with vanishing /* necessarily leads to S = j so that ([3]) is 
eventually approximated by equation (I16p for the streaming particle component, in which one needs 
to have S = j. Physically, one can argue that for large mean free paths A, the streaming particle 
component is supposed to dominate the trapped particle component and all emitted neutrinos 
should directly become streaming. The case S > j would result in even faster decay of /* and an 
effective additional unphysical source of streaming particles in (116p . Therefore, the coupling term 
S in (j23p is limited from above by j. 

Consequently, in the free streaming regime S = j the equation for the trapped particles reads 



df ^ Idlnpjf _ 
cdt 3 cdt dbj 



while the equation for the streaming particles is given by 

(?//''"*■) =^ -§//'*■ 

The two equations are obviously uncoupled. Finally, we note that one would expect the free 
streaming limit to occur outside the neutrinosphere, i.e., with r > Ry m. p7p . 



1.2.7 Reaction regime: S = 

One can argue that the lower limiter E = in (j23p represents the regime in which the roles of 
/* and are changed compared to the streaming regime, reflected by the right hand sides of the 
corresponding equations. Now the trapped particles dominate the streaming particles which can 
only be absorbed. 

The equation for the trapped neutrinos reads 

df ^ Idlnpjf ^. 
cdt 3 cdt duj 

while the equation for the streaming neutrinos takes the form 
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As in the free streaming regime, these equations are uncoupled. The case S < would result in 
an additional source for trapped neutrinos induced by the scattering of streaming neutrinos that 
would eventually lead to an unrealistically big trapped particle component. Note that with S = 
the stationary state of the trapped neutrino equation is given by /* = j/x which is a Fermi-Dirac 
function, i.e., the right distribution function for the trapped particles in thermal equilibrium, cf. 
O p. 822] and [30l p. 1177]. We call the case S = reaction regime since here the reactions ([8]), 
which are encoded in j and X; drive the time evolution of /* together with the second term on the 
left hand side of (1261) . which describes the redistribution of neutrinos in energy space if the gas is 
compressed or expanded (and which also occurs in the trapped particle equation for the diffusion 
and the free streaming regime). Finally, as the diffusion limit we expect the reaction limit to prevail 
inside the neutrinosphere, i.e., with r < R^, in (I17p . 

1.2.8 Mathematical issues 

The main problem of the IDSA is its treatment of transient regions in between the parts of the 
domain where neither the diffusion, the free streaming nor the reaction limit applies. Then, never- 
theless, one of the above regimes occurs. A discussion of this topic by numerical results comparing 
the IDSA with the full Boltzmann equation in the spherically symmetric case is given in [30]. In 
particular, the transition from the diffusive to the free streaming limit in a "semi-transparent" 
regime across the neutrino sphere is quite sensitive and physically essential. The numerical re- 
sults are satisfactory but exhibit the biggest deviations between the IDSA and the full Boltzmann 
equation in these transition regimes. Similar results have been found in a numerical example for a 
model problem in 

In [30] the diffusion limit is derived by a "leading order" Chapman-Enskog-like expansion of the 
Boltzmann equation ^ in which some higher order terms are neglected. However, the solvability 
hypothesis of the Chapman-Enskog expansion ()40p that requires the zero-th angular moment of / 
to be equal to the zero-th angular moment of /o might be violated by physical reasons since /o only 
depends on stationary quantities whereas / represents the conservative evolution of neutrinos in 
time. Further details concerning these issues, e.g., a weakened form of the solvability hypothesis, 
are given in Subsection 12.11 In Subsection 12.21 we present a derivation of the diffusion limit by a 
Hilbert expansion of the time-scaled Boltzmann equation ([3]) without the need to neglect terms. 

The limiters for the coupling term in ()23p can be physically and mathematically motivated (see 
above and Lemma [2.23p . but there does not seem to be a rigorous explanation of these limitations by 
means of the full Boltzmann equation. In Subsection 12. 2^ we obtain these limiters, which represent 
the free streaming and the reaction limits, respectively, by corresponding Hilbert expansions of ([3]). 

Besides the fact that the density range in a core-collapse supernova comprises several orders of 
magnitude, essential processes described by the coupled problem ([I])-(l3]) occur on different time 
scales that, moreover, depend on where in the star the processes take place. Concretely, the reaction 
and collision time scale given by the interaction rates j and x as well as the opacities (po and (pi 
describes the fastest processes of ([3]) at the center of the star. In the semi-transparent regime it 
becomes slower than the transport time scale represented by the second term in ([3|) as well as the 
term F^. In the outer layers it also becomes slower than the hydrodynamic time scale given by 
the terms involving p and f in (l3|). The hydrodynamic time scale is slowest at low density and 
reaches almost c at the center. In contrast, the diffusion time scale is slowest at the center and 
reaches almost c in the streaming regime. We also use time scaling in the Hilbert expansions for 
the derivations of the reaction, diffusion and free streaming limits in Subsection 12.21 
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2 Asymptotic Behaviour of Radiative Transfer and IDSA 



The main purpose of this section is to present different derivations of the diffusion source S — ^ids 
given in Subsection 11.2.51 with the help of asymptotic analysis. In a first approach, we apply 
a Chapman-Enskog expansion to Boltzmann's equation ([3]). Using weakened Chapman-Enskog 
hypotheses, we give an explanation of how to obtain the results for the diffusion limit in pO| App. A] 
within this framework. As a by-product, we also see how the reaction limit in Subsection 11.2.71 can 
be explained in terms of a Chapman-Enskog expansion. 

In a second approach, we apply Hilbert expansions to ([3]) which for the reaction and the diffusion 
limit exhibit strong correlations with the Chapman-Enskog expansions while at the same time they 
require less assumptions. In this framework, we also introduce variations of different scalings for 
time and for the right hand side of the Boltzmann equation. This allows us to quite naturally 
rederive the diffusion limit, but now up to a "clean" O(e^) order of the scaled Boltzmann equation 
with the scaling parameter e > 0. Furthermore, using a Hilbert expansion just for the time scaled 
Boltzmann equation, we derive the free streaming limit as well as the limiter S = j for the streaming 
regime given in Subsection 11.2.61 

2.1 Chapman-Enskog expansion 

The Chapman-Enskog expansion is a well-established method for solving the Boltzmann equation 
for non-uniform gases and in order to derive the Euler and the Navier-Stokes equations as its 
first and second order approximations. It was found independently by S. Chapman [12l [13] and 
D. Enskog [17] in the 1910's. We apply it to our version ([3]) of the Boltzmann equation for neutrinos. 
Our use of the method is based on the presentation in \19\ pp. 115-122] that focusses on Enskog's 
approach. 

If the gas of neutrinos is in equilibrium, the distribution function is given by a Fermi-Dirac 
distribution, see, e.g., [30l p. 1179]. Concretely, this means that / = j/x and = C{f) = 

in ([3]), compare [30 ', p. 1177]. In case of small deviations from the equilibrium, we expect the 
emission and absorption rates as well as the opacities in the collision integral C(/), i.e., the right 
hand side of ([3]), to remain the dominant terms in ([3]) that drive the gas to equilibrium. (This 
reflects the fact that the corresponding processes occur on the fastest time scale.) Therefore, a 
small parameter e > is introduced that emphasizes the weight of these terms on the right hand 
side of the Boltzmann equation, i.e., we write 

W) = e-' {] + J if)) . (27) 

Here, we denote j = ej as well a.s x = and R = eR, or, in particular, (po = £(j)Q and (pi = ecpi in 
the case (f20]l , for the definition of = ejT" according to ([3]) and ([9]) . We assume the terms with a 
bar to be independent of e which for small e results in a small mean free path A = e{x — 0o + 4'i)~^ 
according to (I2ip . Now we look for a solution of (I27p in the form of a power series expansion with 
respect to e, i.e., 

/ = /o + e/i+e'/2 + ... . (28) 

Enskog postulated that the distribution function and, therefore, all /, functions do not depend 
explicitly on time but only implicitly via their dependence on macroscopic observables (and their 
spatial gradients) obtained by velocity moments of / weighted with the collisional invariants. This 
can be justified both mathematically and physically \19\ p. 116]. Here, we do not analyze the 
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analogue of collisional invariants for our case ([3]) but restrict ourselves to the trivial invariant 
function 1 (which is invariant with respect to any variable). Therefore, we only consider the 
velocity moment 

1 



fdfi with 13 = f3{r,uj,t) 



and impose the assumption 



/(r,^,t^,t) = /(r,/i,a;|/3,f 



(29) 



as our equivalent of Enskog's postulate. We can regard f3{r,t,uj) as a measure for the particle 
density of neutrinos for a certain energy u. Its definition and the expansion of / provides f3i, 
i G No, with A = I fidfi such that we formally obtain the series 



/3 = /3o + e/3i + 8^(32 + 



(30) 



and thus the flux defined by 



' dr ' dr'^ 



dt dt dt dt 



Consequently, taking the chain rule into account, we express the partial time derivative of / by the 
formal series 

5/ _ _dl_ d_d^ d(3i df 



k=0 <^ arfe 



i=0 k=0 dr'' 



Inserting the series ([28]) for / in here leads to the Cauchy product 
df dofo 



dt 



dt 



difo _^ dofi 



dt 



dt 



52/0 _^ difi ^ dof2 



dt 



dt 



dt 



+ 



(31) 



where the operators ^ are defined by 



di g d^ d(3., 



d 



dt 



k=0 ^ drk 



0,1,... , 



with the coefficients ^ from the flux expansion. 



(32) 



With the expansion (|3T]) . (l32]l of the partial time derivative operator applied to /, we can now 
expand the operator on the left hand side of ([3]) and write 



(1) 



\(2) 



with 



1=0 



dj-ifi ^ ^dfi_ ^ dfi ^ p dfi ^ ^ dfi 
cdt cdr dr ^ dfi ^ dio 



(33) 
(34) 



using ([6]). The right hand side of ()27p is formally expanded by linearity of J, i.e., 

Jif) = Jifo) + eJih) + e^J{f2) + .... 
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Collecting all terms of the same order i — 1, i = 0,1,..., in e on both sides of the expanded 
equation (f27|l we get the hierarchy of equations 



= j + J{fo), (35) 
T)^f)(i-i) = j^f^)^ i = l,2,... . (36) 

These represent integral equations for the successive terms fi in the power series expansion of /. As 
in the Chapman-Enskog expansion for non-uniform gases, necessary conditions for the solvability 
of (|35p or (j36p are obtained by multiplication of these equations with the invariant function 1 and 
integration over /i. As a result, by (llOp . we get the compatibility equations 

= j-^ r fodf-i^ (37) 



2 



1 



^l'v{f)(^-''^dfi = -^I'j.d^i, i = l,2,.... (38) 

Remark 2.1 If we do not apply the sophisticated operator expansion (I33p but use 'D{fi) instead 
of 'D{f)^^\ motivated by linearity ofT>, the analogous construction leads to the same hierarchy of 
equations as in a Hilbert expansion with scaled reactions and collisions. See Subsection \2.2\ as well 
as Lemma \2.2\ and Remark \2.S\ for the analogies and the differences between the two expansions. 

As mentioned earlier, except for the geometrical relationship (jl7p . the IDSA does not rely on 
information how / behaves with respect to ^. Therefore, it is enough to consider the angularly 
integrated Boltzmann equation ([3]) only, compare also Remark 11.11 In the following, our aim is to 
compute an approximation of the left hand side of this angularly integrated equation, which is also 
known as the total interaction rate s given by 



fl -I fl °^ -. - ,.i oo -. - ,.i 

'■-\L "<fl^^ 4 E (^'-'^(/>'-") - M Lj>-'''^ -i-ilj"^- 

(39) 

In here, and in what follows, we always use formal series and assume their interchangeability with 
integral operators and partial derivations and interchangeability of the latter operators. Further- 
more, we do not use the notation 0{e^) for the A;-th order, > 1, in a rigid mathematical sense 
but only to express that we consider the corresponding series up to the order £^~^ while neglecting 
the rest of the series that starts with the order e^. Since we want to compute s up to second order, 
we need to handle angular means of P(/)W. This can be difficult because of the sophisticated 
treatment of the partial time derivative. Therefore, we base our considerations on the following 
two lemmas. 



Lemma 2.2 //, apart from Enskog's postulate ([29)1 . we assume the strong Chapman-Enskog hy- 
pothesis 



1 

fod^i = /3, (40) 

-1 
1 

f^d^i = 0, i = l,2,... , (41) 

-1 
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then we have 

-1 1 ri 



2 

i.e., after angular integration, the Chapman-Enskog expansion has the same form as the Hilbert 
expansion. 

Proof. By the form of V{fi), cf. i^, and V{fY^^ in (j34p . we only need to prove 



i/;(w)'"-w,)).,.4/;(E^-i)*-o. 



J=0 

i.e., after elimination of the constant c and formal interchange of integration and differentiation, 

1=0 

By ([32]) we obtain ^ = since /3 and its spatial derivatives are assumed to be independent 
variables in our version (I29p of the Enskog postulate. Therefore, all summands in (I32p for A; > 
vanish. Then, hypothesis (j40p provides the statement for i = 0. 

For i > and < / < i we infer by (I32p that 



Here, as for i = 0, all other summands in (I32p for A; > are zero since /3 and, therefore, all /3j do 
not depend on the spatial derivatives of /? by postulate (p9]) . In (|13]) at least one of the terms /3j_/ 
and /3; vanishes by hypothesis (I4ip so that the whole expression becomes zero. Therefore, ()42p is 
trivially satisfied. ■ 



Remark 2.3 The question of unique solvability of the integral equations in the hierarchy (|35p 
and (j36p is the origin of the strong Chapman-Enskog hypothesis (|40p . (j4ip in gas dynamics. Con- 
cretely, in case of non-uniform gases, the compatibility conditions (I37p and (I38p are not sufficient 
for unique solvability of the corresponding original hierarchy equations ()35p and ()36p . Instead, one 
is free to require all the macroscopic information (3 to be already contained in /o, i.e., the strong 
Chapman-Enskog hypothesis, and thus one obtains unique solvability of ()35p and ()36p . /n contrast, 
for the Hilbert expansion, the corresponding compatibility conditions are indeed necessary and suf- 
ficient conditions for the unique solvability of the corresponding hierarchy equations. Existence of 
solutions is provided within the theory of Fredholm integral equations by the compatibility conditions 
that can be regarded as orthogonality conditions for the collisional invariants. Hilbert then proved 
uniqueness of the solutions if initial values are assigned to all Pi, i € Nq. More details can be found 
in 'iJ9, pp. 111-113] for the Hilbert theory and in \1S\ pp. 118/9] for the Chapman-Enskog theory. 
Lemma [KE also holds if 13 is a vector corresponding to collisional invariants since then dS]), i.e., 
the term for k = in the series ()32p , has the form 
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and thus vanishes componentwise, too. We refer to U9\. p. 117] for more details. 
Finally, even if the form of the Chapman-Enskog and the Hilbert expansions are the same, Enskog 's 
hypothesis ([29|) . which is not used in Hilbert expansions, makes a big difference inasmuch each fi 
in the expansion (f28]l of f depends implicitly on every other fi by its dependence on f3, i.e., on all 
Pi, i = 0,1,..., in ([29]). 



Lemma 2.4 We consider the decomposition f = f^ + f^ inducing (3 = + with the angular 
means 

/3* = 5 /_i f^d^i and = \ j\ f'dfx. We assume Enskog 's postulate ([291) and the weakened 
Chapman-Enskog hypothesis 

1 

- J Jodii = i.e., Po = P\ 

°o /I />! \ »D 

i=i -^-1 ^ i=i 

Then with the assumptions 

% = 0{e') and ^ = O^s^) (44) 



we obtain the relations 



\ J'^ Vifpdf, = \ J'^ V{fi)df, + Oie') , iG{0, 1} 



(45) 



and with the assumption 
we obtain 

\ 1' P(/)(°)d/. = \ 1' P(/o)ci^ + 0(e) . (47) 

Proof. As in the proof of Lemma by the form of the time derivative in T>{f)^''^ in 
for (|15]) . we need to prove 

i.e., after ehmination of the e-independent constant c, 

1=0 

for i = 0, 1, which means = ^ + 0{e^) for i = and ^ + %^ = ^ + ©(e^) for i = I. 



By (|32|) we calculate for z = 

doPo _ d/3o d/3o _dPofdl_ dp^\ ^ dPo^^ ^ ^^^2)^ = ^ + o{e'') 



dt dt dp dt \dp dp J dt ' ' dt 
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using the assumptions ()44p on Pq = (3^ and /3* as well as the fact that /3o does not depend on the 
spatial derivatives of /3 in Enskog's postulate ()29p . With the same arguments, also for i = 1, we 
conclude 



5i/3o ^ ^ ^ ^ a/3o 9/3i dPi [3(3 9/3^ 



5t 



dt 



dt dl3 dt dp dt \dp dp 



dPi 



as required. 

The derivation of ()17|) by (jl6|) is the same as the derivation of ()15|) by 
replace ©(e^) above by 0{e). 



for i = if we 



Using the two preceding lemmas, we can write the second order approximation of s as 



1 

s = - 
2 



I / oo 



(i-i) 



.(i-i) 



(i-i) 



\ l'v{fo + efi)df, + 0{e^). 



With this approximation in l\39\i . the equation we want to solve becomes 



V{fo + efi)dfi = ^- 
1 £ 



1=0 ^ -^"1 ^ 



(48) 



which is the sum of the first three equations of the hierarchy (|37p . (j38p in angularly integrated 



form, multiplied by e 



(i-i) 



0, 1, 2, respectively, up to the order 0{e ). 



For further treatment of the left hand side of this equation, the following decomposition of the 
operator V will be helpful. 



Definition 2.5 With the notation as in Subsection \l.l\ we define the operators 



and 



din p 3v 
H 

cdt cr 



dfj, 



dr 



^— + -(l-//2 



'dn 



V 



dlnp ^3v 
cdt cr I cr 



CO 



doj 



(50) 



These operators decompose the operator V = 2?"*" + 'D~ into its symmetric and antisymmetric 
part with respect to ^. In addition, contains all time-dependent coefficient functions. This will 
be used in the next subsection in Theorem 12.191 and Theorem 12.211 

The following result is the fundamental lemma of this and the next subsection. 

Lemma 2.6 Let C{f) be given by (f20l) with x 7^ 0, x + 00 / as well as x + ^Po — 4'i ^- Then 
the first hierarchy equation ()35p has the unique solution 

f - ^ 
Jo — — , 

X 
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which is isotropic and leads to 

V{fT^ = V{fo) (51) 



in the second equation (j36p for i = 1 in the hierarchy, making this equation formally the same as 
in the Hilhert expansion. For the solution of this equation we obtain 

= (^(-^o) + t ^ifo)dl^ + S/^'^i^^ /' ^(/o)/"d/^l • (52) 

X + 00 V X 2 y_i 2 J_i J 



Furthermore, we have 



, V(fo)dfi = - / V+(fo)dfi = ^ + T^u}^ (53) 

2j_i ^■'^^ ^ 2y_i ^■'^^ ^ cdt 3 cdt du ^ ' 



and 

m which = denotes the identity up to D+(D+(-)) terms if we use the decomposition ofD = D^+V^ 
in (Ell) and (ET 



Remark 2.7 The expressions for /q and efi are obtained by taking the zeroth as well as the first 
angular moment of the first two equations of the hierarchy ()35p and ()36p . respectively, and inserting 
the obtained expressions for ^ f^-^ fidfi and ^ f^-^ fipdfi back into (|35]) for i = and (j36l) /or i = 1, 
respectively. Here one uses that the particular form of the truncated Legendre expansion of the 
collision integral C{f) given in equation ()20p is given in terms of f and the zeroth and first angular 
moment of f. The calculations are similar as performed in !30f App. A] and in /5, Sec. 2.4]- 

The arguments of how to obtain (|55p are also sketched in 13(A App. A] and are given in more 
detail in /5, Sec. 2.4]- However, for proper readability of this paper, we need to give a comprehensive 
presentation of the proof of this fundamental lemma here, too. 

In this Lemma, and in the Theorems on the diffusion and the reaction limits to come, /o = j/x 
turns out to be isotropic, so that we obtain the identity 



1 



2 



1 



in which, as in ()15p . the first equation is a slight abuse of notation because its left hand side is 
constant with respect to fi whereas, a priori by integration, the right hand side does not even contain 
a fj,- dependency. 

Proof of Lemma 12.61 Using (jlOp and x ^ 0, the zeroth moment of the first equation (j35p in 
the hierarchy reads 

1 1 /"^ j 

= J - / fodp , i.e., - / fodfi = - . (56) 

2 J-i 2 y_i X 
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Using ()20p . the first moment of ([35]) can be rewritten as 

1 

^ = 2 y 1 -^^^"^^ 

if X + '/'o — 7^ 0- Together with ([56]) and (pOj) in ([35]) . this result leads to 

1 /"^ 1 /"^ j 
= j- xfo- 4>ofo + '/'ot: / fodfJ- + 3/i0i- / /o/xd/i = j - {x + Mfo + (^o- 

2 7-1 2 y_i X 



i-e., /o = j7x if X / and x + (/'o 7^ 0. 

In order to derive (|52p we first prove (|5ip for which, by p4p . we have to show 



dt dt 







(57) 



This property follows from isotropy of /o = j/x which holds since both j and x do not depend 
on fi. Concretely, since /o is isotropic, we immediately have /o = /3o, and since /3 and, therefore, /3o, 
too, do not depend on the spatial derivatives of /3 by postulate (I2HD, we obtain ([^7P from P^ . 

Using ()5ip and (jlOp . the zeroth moment of the second equation (j36p for i = 1 in the hierarchy 
reads ^ ^ 

^ J ^ V{fo)dfi = -x\j ^ efidfi . (58) 

Again, by (|5ip and the special form of the collision integral (j20p . the first moment of (j36p for i = 1 
can be written as 

1 1 3 1 

1^{fo)l^diJ, = -{x + 4>o)^J £fif^dn + (t)i- J i?d^-j efindn 

^ \ J ^^(fo)f^dfi = -{x + (t)o-(t>i)\j ^ehfidfi. (59) 



Now, inserting the expressions (|58p and (j59p for the zeroth and first moment of /i in terms of the 
zeroth and first moment of ^^(/o) into (j36p for i = 1 with the collision integral given by (|2Up . we 
obtain 

V{fo) = -ix + Mefi - t]- f P(/o)rf/i- , \ f V{fo)fidf,. (60) 

X 2 y_i X + '/'o - (Ai 2 y_i 

By the definition (j2ip of the mean free path A, this equation is equivalent to ()52p . 

Now we turn to the proof of ()53|) and ()54j) . Using the decomposition P = V'^ + the first 
equality in ()53p follows immediately from ^ J_iT^ {fo)dfJ. = since 2?^ is antisymmetric in fi and 
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/o does not depend on /u. In order to prove the second equality in (|53|) we calculate 



cdt 



cdt 
df^o 



+ 



dlnp 3v\ 
cdt cr J 



d In p 3,v 
H 

cdt cr 



V 

cr 



UJ- 



2.dfo 

^^^^ 

duj 



dlnp 3v\ V 
cdt cr J cr 



OJ 



dfo 



dp 



1 d In p d(3o 
cdt 3 cdt doj 

Here, the second equality is obtained by formally interchanging integration and differentiation, 
considering /?o = ^ fodp in the first and third summand in the integral, using | f^^ p^dp = ^ 
and the isotropy of /o for the latter and observing that the second summand vanishes completely 
since /o does not depend on p. 

For the approximation "=" in ()54p we also use the decomposition T> = + T>^ and first 
consider the part which, by inserting (152[) . reads 



' V+{Eh)dp = \ 



-1 



+ 3MA^ (P+(/o) + V-{fo))pdp 



\dp. 



Since /o does not depend on p, the term 'D~^{fQ)p in the second inner integral of this expression 
is an antisymmetric polynomial in p and thus vanishes after integration. By the same reasoning, 
the term T)~{fQ)p in the second inner integral is a symmetric polynomial in p so that angular 
integration gives an expression that no longer depends on p. Consequently, the last summand in 
the brackets is linear in p so that the application of the operator (•)) to this expression is 

an antisymmetric polynomial in p that vanishes after the outer integration. By the same reasoning, 
the second summand in the brackets vanishes after the application of this operator and the outer 
integration. With these arguments and ^ 'D^{fo)dp = as seen above we obtain the simplified 
equation 



' V+{eh)dp = l 
-1 ^ 



X + (Po 



X 2 J -I 



V+{h)dp 



\dp, 



(61) 



in which the right hand side only contains terms that undergo the application of the operator T>~^ 
twice. Consequently, this part contains all terms that are neglected in the approximation 

"=" in IPI) and (f5S|). 

Second, the P~ part in the approximation "=" in reads 



V {efi)dp 



-1 



-1 



P+(/o) + P-(/o) + 



^1 



(P+(/o) + P-(/o))d/i 



3/i</>iA- 



{V+{fo) + V-{fo))pdp 



\dp. 
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As for the part, the term 'D^(fo)^ vanishes after integration in the second inner integral. Since 
/o is independent of the term I?^(/o) in the first inner integral is a symmetric polynomial in /i. 
Therefore, since ^ 

/_i T>-{fo)di^ = , the first inner integral does no longer depend on fj, so that the 
application of the operator ( to it is linear in and vanishes after the outer integration. 
The same holds for the first term 'D~^{fo) in the brackets which is a symmetric polynomial in /x so 
that the application of P~ ( (•)) to it gives an antisymmetric polynomial in that vanishes by 
the outer integration. Consequently, we obtain the equation 



V-iefi)df, 











( 


V 




\X + (f>o 





that we can further simplify by observing 
which leads to 



p-(/o) + 3MA 



\dli, 



dr 



(62) 



2^/0 , 



iiAM^ = 0iAP-(/o), 



so that, with (1211). we can conclude 



-1 



{efi)dfi = ^ 



-1 



-1 



:i+0iA)p-(/o 



\df^=l 



-1 



V- {-XV-{fo))dfi. 



X + (Po 

In order to prove the last equality in ()54|] we take (I62p into account and calculate 



i I'v- (AP-(/o)) d/. = i (a/.|^) d^^ 



d 



dr 



dr 



r 



dr 



dn 



1 d / a/3o\ 1 1 2,.dfo, 
3d-r[^^)^2j_rr^'-^^^^''^ 



ld_ 

3 c?r 

]_d_ 



dr 
dr 



+ --X^ 
3 r 3r 



3 9r 



For the third and fourth equality, we used again ^ /^^ fi^dfi = i, the isotropy of /o and A as well 
as /3o = ^ j\ fodfi. The last equality follows from the product rule. 



To finish the proof we note that ()55p follows immediately from the sum of (j53p and (|5 
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Remark 2.8 The approximation (|55p of s in the sense of Lemma \2.6\ is called "leading order 
approximation" in ]3C^ since the "leading ^ 'D^{fo)diJ, term" is included on the right hand side. 
Note that if we introduce an additional scaling for the time variable as in Theorem \2.19l then the 
T)~^(T>^{-)^ terms ()6ip are naturally of order so that (jSSp is replaced by the sum of (153p and 
(j54p in which the first "=" -equation is skipped. This sum is then an approximation of s up to the 
order 0{e^). 

We can now collect the results and formulate our main approximation theorems. 

Theorem 2.9 Let C{f) be given by (pOj) with x / 0, x + '?^0 7^ as well as x + — 01 7^ 0. We 
assume Enskog's postulate (f29l) . Then, with the strong Chapman-Enskog hypothesis in LemmalKM 
the diffusion-reaction-type equation 

dp ^ I dlnp ^ap ]__^L2^^\_j 
cdt 3 cdt duj r"^ dr \ 3 dr J 

is an approximation of the angular mean of Boltzmann's equation (I27p of order 0{e'^) up to 
P"*" (!?"'"(•)) terms in the sense of Lemma \2.6l 

Theorem 2.10 (S in the diffusion limit I) Let C{f) be given by (pOj) with X 7^ 0; X + </'o 7^ 
and X 4'0 ~ (ki 7^ 0. We consider the decomposition f = f^ + f^ inducing (3 = (3^ + 13^ with 
the angular means /3* = ^ f^-^ f^dfi and /3* = ^ J^-^ f^dfi, and we assume Enskog's postulate ()29p . 
Then, with the weakened Chapman-Enskog hypothesis and with the asymptotic properties ()44p of 
(3^ and f3i in Lemma\2.4[ the diffusion-reaction-type equation 



cdt 3 cdt did r"^ dr \ 3 dr J 

is an approximation of the angular mean of Boltzmann's equation (j27p of order 0{e'^) up to 
'D^(T>^{-)^ terms in the sense of Lemma \2.b\ In the same sense 

^..--i^fr^^^Vxr (64) 



dr \ 3 5' 

is an approximation of the (angular mean of the) diffusion source S m ([1 

Proof of Theorems 12.91 and 12.101 First we recall that the asymptotic relation (f48ll is a 
consequence of both Lemma [2.21 and relation (j45p in Lemma [2.41 which follows from ()44p . We insert 
the approximation (j55p into (j48p . Then the equation in Theorem 12.91 follows from /3q = /3. For 
Theorem 12.101 we derive equation (I63p from f3o = (3^ and (3 = f3^ + /3^. The second statement of 
Theorem 12.101 follows from the comparison of (|63p with (jl4p taking into account (|15p . ■ 

In summary, Theorem 12.101 states that the diffusion source (j22p in the IDSA [30, App. A] 
follows from a Chapman-Enskog expansion of second order with the weakened Chapman-Enskog 
hypothesis and the asymptotic properties (j44p of /3* and f3i in Lemma [2. 41 up to T>^(T)~^{-)^ terms 
in the sense of Lemma 12.61 
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Remark 2.11 There is a slight but important subtlety to the conditions imposed on the decompo- 
sition f = + in Theorem \2.10\ and in Lemma \2.4\ In contrast to the same decomposition 
in Section we do not assume /* to he isotropic. Instead, we only impose conditions on the 
angular means /3* and , which reflects the fact that we only make assertions on the angularly 
integrated Boltzmann equation anyway. In particular, we set f3^ = /3o, and we know /3o = j/x from 
Lemma \2.6l Consequently, if /* is isotropic, it is necessarily identical to the equilibrium function 
/o = j/x- This is mathematically not forbidden but physically unrealistic and only satisfied in 
special cases, e.g., in the center of a supernova after long times. Therefore, the assumption in the 
IDSA that /* should be isotropic, i.e, (jl5p . seems inappropriate. As Theorem \2.10\ shows, we do not 
need it in order to establish a diffusion-reaction-type equation as an approximation of the angular 
mean of Boltzmann's equation in the diffusion limit and in order to define the same diffusion source 
as in the IDSA. 

On the other hand, the trapped particle equation ()14p cannot be derived from (|12p if isotropy 
of /* is not assumed since then additional terms occur on the left hand side coming from the non- 
isotropic part of /* . Indeed, in Theorem \2.10\ it is the diffusion term in (|64]l that actually accounts 
for the non-isotropic part of f ^ fo + e/i, which is efi in leading order. If efi is isotropic, we 
have ^T)~ {efi)dfi = since T>^ is antisymmetric in fj, which, by ()54p . entails a vanishing diffusion 
term. 

We close this section by stating the corresponding results for the approximation of the total 
interaction rate s up to first order only, i.e., if we consider the first two hierarchy equations (j35p 
and (|c{6p for i = 1 instead of the first three equations as done above for Theorems 12.91 and I2.1UI It 
turns out that in this case we obtain a reaction-type equation instead of a diffusion-reaction-type 
equation. 

Theorem 2.12 Let C{f) be given by ([20]) with x 0, X + (f>o ^ and x + (po - (t>i 7^ 0- We 
assume Enskog's postulate (j29p . Then, with the strong Chapman-Enskog hypothesis in Lemma \2.2X 
the reaction-type equation 

d/3 ^ldlnp^d/3 _. 
cdt 3 cdt doj 

is an approximation of the angular mean of Boltzmann's equation ()27p of order 

Theorem 2.13 (S in the reaction limit I) Let C{f) be given by (j20P with X 7^ 0; X + 00 7^ 

and X + ipo — 4>i 7^ 0. We consider the decomposition f = f^ + f^ inducing (3 = (3^ + 13^ with 
the angular means /3* = ^ f^-^ f^dfi and (3^ = f^^ f^d^, and we assume Enskog's postulate (I29p . 
Then, with the weakened Chapman-Enskog hypothesis and with the asymptotic properties (|46p of 
13" and (3i in Lemma \2^[ the reaction-type equation 



is an approximation of the angular mean of Boltzmann's equation (j27p of order 0(e). Accordingly, 

^reac •= Xf^^ 

is an approximation of order 0(e) of the (angular mean of the) source term S in ()14p . If we 
additionally suppose f3i = 0{e), then we have (3^ = 0{e^) and the reaction limit proposed in 



Subsection \1.2. 7 



0. 
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is an approximation of order 0{e) ofT, in (jl4p . 

Proof of Theorems 12.121 and 12.131 We first note that the asymptotic relation 

is a consequence of both Lemma 12.21 and relation (I47|) in Lemma 12.41 which follows from (|46p . We 
insert the approximation (j53p into (|66p . Then the equation in Theorem 12.121 follows from /3o = 
For Theorem 12.131 we derive equation (I65p from /3o = /3* and /3 = f3^ + /3^. The second statement 
of Theorem 12.131 follows from the comparison of (j65p with pip taking into account (jlSp . The last 
assertion follows immediately from = /3 — /3* = e/3i + + • • ., the hypothesis x = 0{e~^) and 
the same comparison of equations. ■ 



Summarizing, Theorem 12.131 states that the lower bound for the diffusion source S in the 
IDSA |30l App. A] follows from a Chapman-Enskog expansion of first order with the weakened 
Chapman-Enskog hypothesis and the additional asymptotic properties (|i6]) of and /3i as stated 
in Lemma 12.41 and in Theorem 12.131 Without the latter assumption on /3i we only have = (D{e) 
and the source term Sj-eac is of order 0(1) due to x = However, since on a physical level 

the free streaming component /3* should become very small in a reaction regime, this additional 
assumption might be justified. 

Finally, Remark 1 2 . 1 1 1 also applies to Theorem [2TT3] since here, the angular mean of the Boltzmann 
equation is compared to the trapped particle equation (HH), too. In the reaction limit, S does not 
contain the diffusion term that represents non-isotropic parts of / as indicated in Remark 12.111 
This reflects the fact that in the reaction limit it is physically more realistic for /* to be at least 
close to the isotropic equilibrium function j/x- 



2.2 Hilbert expansion 

In the preceding section, we derived the diffusion limit (Subsections ll.2.5] and ll.2!7|) and the reaction 
limit ( Subsection 1 1 . 2 . 7|) of Boltzmann's equation by means of Chapman-Enskog expansions. In this 
section, we will use Hilbert expansions, first in order to derive the diffusion and the reaction 
limits with less assumptions on the asymptotic behavior of the components in the expansion. 
In fact, it turns out that these additional assumptions (as given in Lemmas 12.21 or [2^ respectively) 
were only needed above to obtain a Hilbert expansion setting from a Chapman-Enskog expansion 
ansatz. Moreover, in a second step, we can even strengthen an assertion by quite naturally deriving 
approximations of a certain e-order for the diffusion limit without the restrictive exception of 
"up to P+ (!>"'"(•)) terms" as in Lemma 12.61 and Theorems I2.9f[2.131 This will be achieved by an 
additional scaling of the time variable in the operator V as done, e.g., in Degond and Jin jl5j . 
compare also [1], [20] and [39] and the literature cited therein. Finally, by just considering the 
time scaling without the scaling of reactions and collisions, we can also derive the free streaming 
limit of the IDSA introduced in Subsection 11.2.61 For an introduction to Hilbert expansions we 
refer to [191 pp. 109-115]. The ansatz is the same as for Chapman-Enskog expansions since it starts 
with a scaling of the right hand side of Boltzmann's equation as in (I27p and an expansion of the 
distribution function / in powers of e as in (128p . However, instead of the sophisticated operator 
expansion (133p that is based on the Enskog postulate (I29p . Hilbert uses linearity and the formal 
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continuity of T> which leads to the simple expansion 

Vif) = V{fo) + eV{h) + e'Vih) + .... (67) 

In addition, instead of the Chapman-Enskog hypothesis (j40p and (j4ip . Hilbert assumes another 
condition that leads to unique solvability of the obtained hierarchy of equations which shall not be 
discussed here. (With regard to this topic we refer the reader to 1191 p. 111].) 

Using the expansions (128p for / and (I67p for T> in Boltzmann's equation ()27p we obtain the 
hierarchy of equations 

= j + J(/o), (68) 
V{fi-i) = J{fi), i = l,2,... , (69) 



corresponding to e-powers i — 1 for i = 0, 1, ... . With the scaling in (I27p that is applied only to 
the reactions and the collisions, we now derive Theorems I2.91[27T3] and, in particular, the diffusion 
and the reaction limit of Boltzmann's equation in the same way as based on Chapman-Enskog 
expansions. Remark 12.111 also applies here. 



Theorem 2.14 Theorem \2.9\ also holds without assuming Enskog's postulate (1291) and without the 
strong Chapman-Enskog hypothesis in Lemma \2.S\ 



Theorem 2.15 (S in the diffusion Umit II) Theorem \2.1(A also holds without assuming En- 
skog's postulate (|29p . without the weakened Chapman-Enskog hypothesis in Lemma \2.4\ and without 
assuming the asymptotic properties (I44p of /3* and (3i . 



Theorem 2.16 Theorem \2. 12\ also holds without assuming Enskog's postulate (|29p and without the 
strong Chapman-Enskog hypothesis in Lemma \2.2i 

Theorem 2.17 (S in the reaction limit II) Theorem \2.13\ also holds without assuming Enskog's 
postulate , without the weakened Chapman-Enskog hypothesis in Lemma \2.4\ and without as- 
suming the asymptotic property (jl6|) of f3^ . 

The proofs of Theorems l2.14fE7T71 can now be based on the hierarchy of equations (|68p and (|69p 
according to Hilbert instead of the hierarchy of equations (|35p and (j36p according to Chapman and 
Enskog. Lemmas (|2.2p and (j2.4p were only established to be able to use the former if one starts 
with the latter. Therefore, the assumptions in these lemmas are not needed here, and the other 
parts of the proofs are exactly the same as above since they are anyway based on the operator 
expansion (j67p . These arguments make it evident that for our purposes, Hilbert expansions seem 
to provide a more natural approach than Chapman-Enskog expansions. 

In the next result, we get rid of yet another restriction, now in the assertion of Theorem 12.101 
or 12.151 in which the approximations are in fact of order O(e^) without the need to neglect any 
further terms. This result will quite naturally be obtained by the introduction of a time scaling 

t = (70) 

in addition to the scaling of reactions and collisions in ()27p . compare [15]. Concretely, we consider 
the scaled Boltzmann equation 

V,{f)=e-'{j + J{f)) (71) 
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where represents the operator D with scaled time (j70p . As in (j27p . we have 

. i ~ X o ^ 

j = -, x = - , R = — , 

e e e 

in which the latter two quantities give J = J je according to ([3]) and ([9]) and in which R = R/ e \s 
induced by 0o = <^o/^ and (/>i = c^i/e in the approximate case ([20|) . Again, we assume the terms 
with a bar to be independent of e. 

By the time scaling ([70]) . the velocity v = dx/dt of the background matter is scaled as 

V = ev 

which by ()49p and ()50p results in a scaling of the operator T) as 

= eV+ + V~ (72) 

since the time dependency is completely contained in the denominator of the coefficients in the 
part T>^ of the operator T> which is symmetric with respect to 

Remark 2.18 By this time scaling, we consider big reaction and collision rates over long times and 
small velocities of the background matter, which results in a smaller influence of the aberration and 
the Doppler effect arising from the motion of the matter described by the part P"*" of the Boltzmann 
operator D. We do not discuss here whether this additional scaling of time is physically reasonable 
for supernova neutrinos in the diffusion regime for which we consider it. In order to decide this, 
we would need a detailed scale analysis of the Boltzmann equation for radiative transfer that should 
be performed elsewhere. However, the same scaling has already been successfully applied for the 
Boltzmann equation in gas dynamics J 15^ - We also refer to 14, Rem. 5] where for the diffusion 
approximation of a transport equation with the same scaling structure as ()7ip with (I72p . the scaling 
of the mean free path and of the time has been rigorously justified. Analogously scaled kinetic 
equations are, e.g., considered in 120^ and ]3iA Sec. 5.1] with respect to the diffusion limit (without 
the time scaling one rather obtains the hyperbolic limit, see, e.g., /gPl Sec. 5.1] or \W]). Finally, 
also in our case of radiative transfer, the additional time scaling turns out to be well suited for the 
derivation of the diffusion source as we will see in the following. 

Inserting the expansion ()28p for / into ()67p with the scaled operator V = leads to the formal 
series 

-Deif) = V-{fo) + e(P+(/o) + + e^V+ifi) + ^-(/s)) + . . . . (73) 

With this series as the left hand side of the time and reaction scaled Boltzmann equation ()7ip , 
we obtain the hierarchy of equations 

= j + J{fo), (74) 
-D'ifo) = Jifi), (75) 
V+{f,^2) + V-{f,^i) = Jifi), i = 2,3,..., (76) 

that corresponds to e-powers z — 1 for i = 0, 1, . . . . 

The fact that the components in expansion (I73p corresponding to the operator part P"*" appear 
in a higher order than the ones corresponding to V~ allows us to derive the diffusion-reaction-type 
approximation of the Boltzmann equation as well as diffusion source up to the order 0{£^) without 
having to neglect D+(D^"(-)) terms since they are now naturally of higher order. 
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Theorem 2.19 (S in the diffusion Umit III) Let C{f) be given by ([20]) with x 0, x + 4>o ^ 
and X + </'o ~ 9^1 7^ 0. We consider the decomposition f = + inducing (3 = (3^ + 13^ with the 
angular means /3* = ^ J^^ f^dfi and /3'^ = ^ J^^ f^dfi. Then the diffusion-reaction-type equation ([63]) 
is an approximation of the angular mean of Boltzmann's equation (|71|) of order 0{£^). Furthermore, 
expression (|64p is an approximation of order ©(e^) of the (angular mean of the) diffusion source S 
in (IE 



The proof reuses arguments by which Lemma 12.61 was proved. In fact, it reUes entirely on the 
variant of Lemma 12.61 that is based on the Hilbert hierarchy ()74p - (j76p instead of the Chapman- 
Enskog hierarchy (j35p and (j36p . This variant is proved in the following. 

Proof of Theorem 12.191 The first equation (|74p of the hierarchy (j74P " (j76p is the same as 
the first equation (jSSp in the Chapman-Enskog hierarchy (j35p and (|36p considered in Lemma 12.61 
and, therefore, has the unique solution /o = j /x for which (|53p holds. 

The second equation ([75]) is treated as equation p6|) for i = 1 with ([5T]) in the proof of 
Lemma 12.61 We only need to replace T> by T>~ in ()58p - (|60p and in the result (j52p . The for- 
mulas even simplify since by the first equation in (I53p the right hand side in (]58p . i.e., the angular 
mean of /i, and the second summand in (I52p vanish so that the latter now reads 

efi = (v-{fo) + ■il^<Pi\\ C V-{fo)fidfi) . (77) 

X + (Po \ 2 7_i J 

By the proof of Lemma 12.61 it becomes clear that with this result the assertions (154p and (155p are 
even true if we replace "=" by "=" since P+(X'+(-)) terms no longer occur. 

To finish the proof we add the angular means of the first three equations in the hierarchy 
(j7i|) - ([76|) and, considering again (fTOj) . obtain 

I |\2?(/o) + 2?-(e/i))d/x = i + | f^(j2e'f}jdfi 

taking the formal limit in the series and considering that j and x are of order 0{e~^). The required 
diffusion-reaction- type equation (f63l) as an approximation of (ITT]) of order O(e^) is obtained by 
inserting ()53p and ()54p on the left hand side while setting /3o = /?* and /3 = f3* + (3^. The corre- 
sponding approximation of the diffusion source (164p follows from the comparison of (163p with (jl4p 
while taking into account (jlSp . ■ 



For more details concerning the statement on the diffusion source in Theorem 12.191 compared to 
its appearance in the trapped particle equation [T3] we refer to Remark 1 2 . 1 1 1 that also applies here. 

Remark 2.20 If one treats equation ()76p for the order , i.e., 

P+(/i)+p-(/2) = J(/2) (78) 

according to (j58p - ()60p . one obtains an analogous formula as (j52p for e/2. Concretely, one only 
needs to replace I^(/o) by ^^(/o) -\-'D~(efi) on the right hand side of (j52p . Plugging this into (j78p 
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and taking into account that fi in (j77p is only made up by T>~ terms, one can see that there are still 
no terms as in (I6ip occurring in the second order equation of the hierarchy (I74p - ()76p . 

These terms only appear in the equation for the order since this one contains the term D~^{f2) 
and the formula for e/2 is also made up by a P+(/o) term as just seen. 

The reaction-type equation (j65p is not available with the time scaling (j70p since the angular 
means of the first two equations (j74l) and (j75p only provide the zero operator on the left hand side. 
In order to achieve this with the first two equations of the hierarchy, one would need the (inverse) 
operator scaling T)^ + eT)~ instead of (172p . This could be achieved with the inverse time scaling 
t = et (i.e., short times considered) while scaling the right hand side by e^^ (i.e., very large reaction 
and collision rates considered). However, the third equation (j76p of the hierarchy would then be 
made up entirely by terms on its left hand side. 

Conversely, if we consider very long times, i.e., the time scaling t = tje^ with the reaction 
and collision scaling as in ()27p . then the first T>^ term, in particular the time derivative, only 
appears in the fourth equation of the hierarchy. As in the proof of Theorem \2.19l we then obtain 
only the diffusion term in (I54p by the first three equations, i.e., a stationary diffusion equation as 
an approximation of the corresponding scaled Boltzmann equation of order 0{e^). However, if we 
compare this equation with (I14p . we obtain the reaction term in the source term S. 

Finally, there is also a version of Theorem \2.19\ based on a Chapman-Enskog expansion and 
Enskog's postulate ()29p instead of a Hilbert expansion. In this context, the time derivative in the 
term P"'"(/j_2) in (j76p would be replaced by the corresponding order term as in (|3ip. To establish 
the result corresponding to Theorem \2.19\ for the Chapman-Enskog expansion with the same time 
scaling ()70p . we need to force that term to be close enough to the time derivative term in D'^(/j_2) 
for i = 2 as done in Lemma \2.4\ Concretely, we need assertion (j47p in which 0{e) is replaced 
by C(e^) in order to relate the angular mean of the operator in the Chapman-Enskog expansion 
with the one in the Hilbert expansion up to the correct order. To guarantee this, the proof of 
Lemma \2.4\ shows that we need assumption (j46l) in which 0{e) is replaced by 0{e'^). 

In order to obtain the free streaming limit we need to consider long times while assuming 
moderate reaction and collision rates. Therefore, we consider the time scaling (j70p . i.e., the operator 
scaling (I72p . and no longer impose the scaling on the reaction and the collision terms, thus we deal 
with the scaled Boltzmann equation 

eV+{f) + V-{f)=j + J{f). (79) 

With this scaling of the Boltzmann equation and the expansion ()73p of the left hand side, 
collecting the same e-powers, we arrive at the hierarchy of equations 

^"(/o) = j + J{fo), (80) 
D+(/i_i) + p-(/i) = Jifi), i = l,2,..., (81) 

collecting the terms of the same e-powers e*, i = 0, 1, 2, ... . The first equation (jSOp of this hier- 
archy is not time-dependent and turns out to provide the stationary state limit for the streaming 
component of /. Since only the angular mean of ()80p is considered, we do not need to assume 
the special form (j20p for the collisions since by (jlOp they do not play any role here. However, in 
contrast to (I62p and the derivation of ()54p . we cannot assume isotropy of /o here. 
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Theorem 2.21 (S in the free streaming hmit) We consider the decomposition f = f'' + 



inducing (3 = + {3^ with the angular means /3* = 2 j_i f^dfi and fi^ = 2 I-i f^d^ as well as the 
expansions (j28p and (j30p for f and its angular mean (3. Then the stationary state equation 

1 g / 1 ri 



2I 

r - 
2 



(82) 



is an approximation of the angular mean of Boltzmann's equation (I79p of order 0(e) with = /o 
and thus /3* = /3o and (3^ = (3 — (3^ = 0{e). In addition, 

Sfree •= J 

is an approximation of order 0{e) of the (angular mean of the) diffusion source S in (jl6p . 



Proof. By (jSOp and (|10p . the angular mean of the first equation (jSOp in the hierarchy ([80 
and (i8T]l gives 
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-1 
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fodjd. 
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Integrating by parts we get 
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fotJ'dfJ, : 
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so that, after interchanging integration and differentiation, (|82p can be written as 



a 1 



2 1 



r 2 



By the product rule, this equation can be more compactly formulated as 



,1 



Now, setting = /o and thus /J'' = /3o provides equation ([82]) and the first statement of Theo- 
rem |2?2I1 The second statement on the source term S follows from the comparison of (|82p with (|16p . 



Remark 2.22 The scaling in ([79p contains the time scaling (1 70 p . i.e., one considers long times, 
and no reaction and collision scaling. On a physical level, one can argue that one should in fact 
introduce the scaling 

j = £j , X = ex , R = eR, J = ej 

of reactions and collisions because in the free streaming limit, e.g., far away from the core of a 
collapsed star, one can assume vanishing reaction and collision rates. Since in the free streaming 
limit we also assume the distribution function to be in a stationary state, we can then consider an 
even larger time scale, e.g., the time scaling 



t 
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The resulting scaled Boltzmann equation then reads 



(83) 



and exhibits the hierarchy 

0, 

3 + J{k) , 

Jifi^i), i = 2,3,..., 

corresponding to e-powers e*, i = 0, 1, . . . . Considering the sum of the angular means of the first two 
equations of this hierarchy and the same arguments as in the proof of Theorem \2.21l one can then 
derive equation (\82h with = /o + ^/i as an order 0{e'^) approximation of the scaled Boltzmann 
equation (j83|) . Analogously, the source term j is then also determined up to the order O(e^). 
In fact, one obtains ©(e") approximations for any n £ N i/ one considers the scaled Boltzmann 
equation ([55]) with in the time scaling replaced by e" and e in the reaction and collision scaling 
replaced by e^~^ . 

Note that we do not require the collision term C{f) to be of any specific form here. Since we 
only take angular means of the equations in the hierarchy, by (jlOp it is even irrelevant how we scale 
the collisions. 

Finally, we remark that one can also consider these scalings for the free streaming limit in 
Chapman-Enskog expansions with Enskog 's postulate (I29p . In fact, since in contrast to the reaction 
and the diffusion limit, the time- dependency is always of higher order in the free streaming limit, 
the Chapman-Enskog expansion is identical to the Hilbert expansion up to the order that needs to 
be considered here. 



The following table summarizes the scaling results in Theorems 12.171 12.191 and 12.211 



~~ — — time scaled 
reac./coll. scaled ~" — 


yes 


no 


yes 


diffusion limit 


reaction limit 


no 


free streaming limit 


full Boltzmann 



V'if^) = 
V+(f,^2) + V-(fi) = 



With the source terms S^iff, Sreac and Sfree in the diffusion, reaction and free streaming limit 
by Theorems 12.171 12.191 and 12. 2H we can now define a global source term 



-'glob • 



mm < max 



mm < max 



■ 1 d 



1 ^free^ 

A a/3*' 



r"^ dr \ 3 dr 
or, with Sreac instead of Sreac, the IDS A source term 



-'glob 



mm < max 



^diff 



free 



mm < max 



]_d_ 



(84) 



(85) 
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as in (p3|) . However, in contrast to the system of coupled equations p6j) and (|23|) of the IDSA, 
we do not need to assume isotropy /* = /?*, i.e., (fTSj) . so that we consider the equations 

^ + i^^M = J- _ _ S (86) 
cdt 3 cdt 5a; 

given by Theorems 12.171 12.191 and 12.211 and coupled by either S = Sgiob or S = Sgiob ■ 

We close this investigation by a short motivation why the choice of limiters from above and 
below for the source term Sdifr as in the min-max-expressions (I84p and (IBSp is indeed mathematically 
reasonable. We base our considerations on the approximation (|86p of the angularly integrated 
Boltzmann equation with as yet undefined but, clearly, isotropic S. Recall that we have X = j + X 
by Subsection 11.1.21 

Lemma 2.23 Let the background matter, i.e., p, the reaction rates j and x o,nd the source term T, 
in ([86]) be in a stationary state limit, i.e., independent of time. Then S is limited from above and 
below by 

-X<^<3- (87) 

Proof. Since p is independent of time, the second term on the left hand side of ()86p vanishes 
and we obtain the reaction equation 

cdt 

with time-independent data j, x S on the right hand side. Consequently, this equation possesses 
the analytical solution 

/?*(t) = /3*(0)e-^** + (1 - e-^**) . 

X 

For t — 7> oo the solution approaches the stationary state limit 

/3^= lim /3*(t) = :^. 

t^oo X 

Since, just as /*, the function /3*(t) = \ f^i f^{t)dp. is a distribution function for every t > 0, 
so is i.e., we have < < 1 and, therefore, 

i - S 

< < 1 . (88) 

X 

With X = J + X this estimation is equivalent to (f87|) . ■ 

By the stationary state assumption in Lemma 12.231 it turns out that the source term Sfree = j 
in the free streaming limit is indeed an upper bound for S because of the non-negativity of the 
distribution function. The bound < S in the IDSA (j85p is more restrictive than — x < S in (|87p . 
and even more so xf^^ < S in (I84p . However, — x < S in (I87p arises from the upper bound /3* < 1 
which one can regard as quite extreme. If instead, one uses the thermal equilibrium function j/x 
(cf. [30l p. 1177]) as an upper bound for i.e., if one substitutes 1 by j/x in (|88]l . one obtains 
< S. 
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